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ABSTRACT 

The polarization potential plays a decisive role in shaping up the cross sections 
in low-energy positron collisions with atoms and molecules. However, its inclusion 
without involving any adjustable parameter, is still a challange. Here we summarise 
various other techniques employed so far for positron collisions and discuss a new, 
nonadjustable and very simple form of the polarization potential for positron-atom 
(molecule) collisions below the threshold of positronium formation. This new po- 
tential, recently proposed by us, is based on the correlation energy (6corr) °f a 
single positron in a homogeneous electron gas. The e corr has been calculated by 
solving the Schrodinger equation of the positron-electron system and fitted to an 
analytical form in various ranges of the density parameter. In the outside region, the 
e C orr is joined smoothly with the correct asymptotic form ( — > w ^ ere a o is the 

polarisibility of the target ). We tested this new positron correlation polarization 
( POOP ) potential on several atomic and molecular targets such as the Ar, CO, 
and CH 4 . The results on the total and differential cross sections on these targets 
are shown here alongwith the experimental data where ever they are available. 


L Introduction 

In the positron(e + ) scattering with multi-electron atom ( or molecule ), a true e + polarization 
potential is very difficult to incorporate without involving any adjustable parameter. Two simple 
stratagems have been quite popular: one, to use the electron polarization potential as such for the 
corresponding positron collisions and two, to employ a phenomenological form under the tuning 
procedure. However, both approaches are unsatisfactory and usually fail at the differential cross 
section (DCS) and annihilation parameter ( Z c // ) level; although one may be succesful in getting 
good agreement for the integral quantities. 

It is only recently that several theoretical attempts have been made to consider the polar- 
ization of the target atom ( or molecule ) by the e + at the ab initio level^ 1 8 ); however, these 
rigorous calculations are not totally parameter-free and suffer from including near-the-target non- 
adiabatic effects and also the correct values in the asymptotic region. Although, the question of 
nonadiabaticity may not as crucial for the positron projectile as it is for the electron case ( see 
later ). The polarization effects dominate the scattering process at very low energies ( below about 
5 eV ). At somewhat higher energies ( roughly above 5 eV ), the DCS s are still quite sensitive to 
such charge distortion effects. Here we are concerned only at low energies, particularly below the 
positronium formation and/or any electronic excitation thresholds. For the one— electron positron- 
hydrogen case, the issue of polarization potential has recently been discussed by Abdel -Raouf 9 . 
A very recent and comprehensive review 10 on the e + -molecule scattering gives details on various 
approaches used so far to include polarization effects in positron scattering. 
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It is, therefore, always desirable to find some kind of simple form of the polarization force 
without involving any fitting parameter. At low impact energies, where the impinging positron is 
moving very slowly relative to the target electron, the distortion of the electronic charge cloud is 
quite different than the corresponding electron impact case. Asymptotically, the polarization po- 
tential can be assumed to be the same for both the projectiles ( which is true up to the second-order 
perturbative theory only ); however, at near— the— target encounters, both the charged projectiles 
interact differently with the charge cloud of the target. In this talk, we discuss a new parameter-free 
form of an approximate e 4 " polarization potential which is simply a function of the target charge 
density and in addition is energy-independent. 

A summary of earlier usage of various models for polarization potential in positron-molecule 
collisions is given by Morrison and coworkers (1_2) . Recently, Elza et al 8 have investigated various 
aspects of polarization/ correlation effects in low-energy positron-N 2 collisions via a two-parameter 
model. Tennyson 3 and Tennyson and Morgan 4 have applied the R-matrix technique to positron- 
molecule (H 2 , N 2 and CO ) collisions. However, the inclusion of polarization force in all these the 
so called ab initio methods is not either satisfactory nor complete. From numerical point of view, 
these procedures are not easy to apply for a general positron— molecule system. The collision of 
positrons with polyatomic targets is even more complicated. We^ 11-12 ^ have reviewed the situation 
on the polyatomic molecules with respective to the polarization effects and the comparison with 
the experimental data. 

In the next section, we describe the new positron polarization potential and in section III, the 
numerical techniques are summarised. In order to demonstrate the success of the new positron 
correlation polarization (PCOP) potential, in section IV, we present some calculations on the 
differential and integral elastic cross sections for the atomic ( Ar ) and molecular ( CO and CH 4 ) 
targets below the positronium ( Ps ) formation threshold ( i.e. E < 10 eV ). Concluding remarks 
are given in the last section. 


II. The New Positron-Correlation-Polarization (PCOP) Potential 


as 


Asymptotically, for a general e + -molecule collision system, the polarization potential behaves 


v pol{r,0,<fi ) = ^-[o:o(47r)^ S® 1 + a 2 (^7r)? S® 1 + S 2 21 ], 


( 1 ) 


where the S™ g is a real spherical harmonic ( see Ref. 13 for its definition and various properties ), 
( r ,9,<f>) are the coordinates of the projectile referring to the center of the target and the spherical 
(oo) and nonspherical (a 2 and a \ ) polarisibilities are expressed in terms of the polarisibility tensor 
an of the target, namely, 

1 2 11 
<*o = ^( a n + <* 22 + <* 33 ); «2 = 2(033 - - 2 ^ 22 ); o 2 = a u - o 22 . 

The above form ( Eq. 1 ) of the polarization potential is accurate at large r values up to the 
second-order perturbation theory. The problem arises when the projectile is near the target. A 
simple way to remedy this difficulty has been to multiply Eq. (1) by a cut-off function depending 
upon some adjustable parameter; however, this approach is unsatisfactory , although the results may 
be forced to agree with observations ( see for example, Darewych 14 , Horbatsch and Darewych 15 ). 
For positron collisions, most of the calculations prior to 1984 used an electron polarization potential 
(EPP) assuming that such polarization effects are not sensitive to the sign of the charge of the 
projectile. Morrison and his group ( Morrison et al 1 ; Morrison 2 ; Elza et al 8 ) strongly advocated 
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for a. need to generate a true positron polarization potential (PPP) rather than employing the 
EPP; they strengthened their point by presenting detailed calculations on the posit ron-Ho and N 2 
systems and comparing them with experimental data. Although earlier calculations on the positron 
collisions using the EPP gave good results as compared to experimental <7^ values, however, these 
theoretical results are generally poor at low energies and not qualitatively good at all for differential 
cross sections at any energy (16-17) . Unfortunately, the more rigorous calculations based on the 
variational polarised— orbital theories are not satisfactory either; Elza et al have to introduce a 
cut-off function in both the short and long ranges and adjust two parameters to make theory and 
measurement in close agreement. Even in a more sophisticated R— matrix approach, an accurate 
inclusion of polarization effects has not been achieved yet^ 3 5 h 

Here, our goal is to look for a computationally simple form of the positron polarization potential 
which is different from the corresponding electron potential and virtually free from any adjustable 
parameter. The basic philosophy of the present approach is similar to the method of 0 Connel 
and Lane 18 for the case of electron scattering based on the correlation energy of the target in the 
presence of an incoming electron. The present positron polarization potential is also based on the 
correlation energy of a localized positron in an electron gas and its hybridization with the correct 
asymptotic form. Here, we think the incoming e+ as a charged impurity at a fixed distance in 
an homogeneous electron— gas. In positron annihilation experiments, a fundamental question to be 
asked is how the electron— positron interaction distorts the electronic structure of the system under 
investigation. 

The e + correlation energy e corr in an electron gas has been evaluated phenomenologically (19-20) 
as well as using the Bethe— Goldstone type approach 2 *. Recently, Arponen and Pajanne have 
applied a completely new approach to the problem of a light impurity in an electron gas. In their 
method 22 the electron gas is described by a set of interacting bosons representing the collective 
excitations of the random— phase— approximation (RPA). Very recently, Boronski and Niemmen 
have described the density functional theory of the electron-positron system and presented the 
results on the positron— electron correlation energy as a function of the density parameter r a ( see 
later ) for different n+(r)/n_(r) ratios including the case of one positron in a homogeneous electron 
gas. Here 77 ,- 1 - and denote the densities of positrons and electrons respectively. 

The physical picture of the positron correlation in an electron gas is as follows. When the 
incoming positron enters the target electronic charge cloud, we can assume the positron as localized 
instantaneously and correlating with the surrounding electrons of a given density, n_(r). The wave 
function of the positrons in such an electron— positron plasma, can be written as 


-L 2 t + ( r) + [/i* c (n+(r)) - ^(r) + l)^(r) - e+tb?(r), 


here f. i xc is the exchange— correlation potential ( which is zero in the present one positron case 
), <£(r) is the Hartree-Coulomb potential and E e c ~ p is the positron-electron correlation energy 
functional. The Eq. (2) has been solved numerically in a self-consistent manner 23 . Based on the 
paper of Arponen and Pajanne 22 , Boronski and Nieminen 23 have given explicit expressions for the 
positron-electron correlation energy, e corr (r s ) interpolating it for the whole radial region. These 
expressions are obtained without giving any divergence problems in the calculations of annihilation 
rates over the entire range of the density parameter r a ( Kallio et al 24 ). In their work, Arponen 
and Pajanne 22 have developed a new approach to solve the problem of a charged impurity in an 
electron gas. The correlation energy, e corr , is calculated from the ground-state expectation value 
of the Hamiltonian which describes the electron gas plus the incoming positron fixed at a distance. 
In the evaluation of e corr , the positron-electron interaction has also been considered ( see Eq. 2 ). 
The analytic interpolated expressions for the e corr in the whole range of the density parameter r 3 
( | 7 rr 3 p(r) = 1, where p( r) is the undistorted electronic density of the target ) are given as follows: 
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2e corr (r s ) = -0.92305 - 0 ; 


t 3 < 0.302, 
0.302 < r s < 0.56, 


. , , 13.15111 2.8655 

2e C orr(r a ) — — - — — + — ; — — — 0.6298; 0.56 < r 3 < 8.0, 


(r* + 2.5) 2 (r, + 2.5) 


and finally, 


(3a) 

( 34 ) 

(3c) 


2e C orr(rc(r 4 )) — 179856. 2768n 2 -f- 186.420(77 — 0.524; 8.0 ^7 * 4 < oc, (3d) 

where 77(r 4 ) is the electronic density corresponding to the density parameter r 4 . 

The PCOP potential, defined as a functional derivative of the correlation energy with respect 
to p(r), can be derived conveniently from the following equation in terms of functional derivative 
of the density parameter 24 , 

Tcorr(r) = ^l--r,— 


Finally, we obtain the following form of the F corr (r) ( in atomic units ) from Eqs. (3)-(4): 
for t, < 0.302, 



2F corr (r) = 


-1.30 

—— + (0.051 ln(r 4 ) - 0.115)ln(r 4 ) + 1.167; 
v r * 


(5a) 


for 0.302 < r, < 0.56, 

2F corr (r) = -0.92305 - ( 56 ) s 

. . 

and for 0.56 < r, < 8.0, 


2F corr (r) = - 


8.7674r 4 
(r, - 2.5) 3 


+ 


-13.151 + 0.9552r 4 
(r 4 + 2.5)2 


2.8655 


(r 4 + 2.5) 


— - 0.6298. 


(5c) 


Note that for molecular systems the short-range e corr (r 4 ) is to be divided by a factor of 
(2 £ + 1)/ \/47r to account for molecular orientation not considered by Arponen and Pajanne 22 . Here 
we do not worry about the 8.0 < r 4 < oo region, as this range is beyond the crossing point where 
the polarization potential is accurately described by the asymptotic term ( Eq. 1 ). It is to be 
noted that the interpolation formulae for the correlation energy ( Eqs. (3) ) were formulated in 
such a way that for the limit r, -too , the e corr reaches the value of Ps _ ion energy, i.e., -0.262 
a.u.. In the present positron scattering case, we realize that in the r, — *• oo limit, the correlation 
energy approaches the correct asymptotic form of the polarization potential ( the same prescription 
as suggested by O’Connel and Lane 18 ). 

Thus, the PCOP interaction potential, F p ^f op (r), for the e + -molecule/atom system is given 

y p F I OP M = V„„(r), r < r C) (6a) 
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and by Eq. (1) for the r > r c range. Here r c is the radius where the V corr and -a 0 /2r 4 ( or 
a 2 /2r' 4 ) terms cross each other for the first time. In addition, we will also report the similar 
cross sections under the ECOP potential, which has recently been employed for positron-molecule 
scattering ( Jain (26 ~ 27) , Gianturco et al 28 ). Even though the EPP results are encouraging for 
some molecular targets ( Jain and Thompson 29 ), we feel that it is more appropriate to find a true 
positron polarization interaction. 

The new PCOP potential ( Eq. 6 ) has several favourable points worth mentioning here: 
first, it involves a true correlation of the incoming positron with the target electrons at short 
distance encounters and exihibits correct behavior in the asymptotic region; second, it is very easy 
to calculate and convenient to incorporate into any model potential approach; third, it is quite 
different from the corresponding EPP and finally, ( see later ), it gives qualitative good results for 
the total cross sections for several atomic and molecular targets as compared with experimental 
data. 


III. Scattering Parameters 

In fact, the numerical techniques to solve the scattering equation for the e~*" wave function 
are standard as employed for the electron scattering case. For the atomic target ( Ar ), we use 
the variable— phase— approach (VPA) 39 in order to determine phase— shifts at each energy, more 
details of the VPA approach are given in Ref. 31. The optical potential of the e+-Ar system is 
determined very accurately from the numerical Hartree-Fock wave functions of the target 32 . In 
order to preserve numerical accuracy, convergence tests were carried out with respect t.o^ radial 
integration and number of partial waves retained in the evaluation of various cross sections 

For the diatomic molecule (CO) case, the scattering equations are set up in the single center 
formalism under the body-fixed (BF) adiabatic-nuclei-approximation (ANA). The final coupled 
scheme is formulated in the integral equation method 33 . More details for the positron-CO calcula- 
tions can be found in our earlier (26_27 ' and recent 17 papers. 

Assuming the CO molecule in its ground electronic ( l<r 2 2<r 2 3<r 2 4a 2 l7r 4 ; 1 E ) and vibrational 
states, the equation of the continuum positron wave function, P(r) in the single center formalism 
under the BF ANA can be written as, 

[V 2 + k 2 - 2V (r)]P(r) - 0, (7) 

where k 2 is the positron energy in Rydbergs and the interaction potential V(r) includes the repulsive 
static and attractive polarization forces. Expanding the E(r) in terms of Legendre projections, v \ , 

ax 

E(r) = E, t (r) + V^r) = [vf T (r) + uf OL (r)]P x (cos ^), (8) 

A=0 


we obtain the following set of coupled differential equations for the continuum function P( r) for a 
given symmetry A, 

l|r - = E M, (9) 

where the potential matrix Vn< is determined as usual 34 . Here A corresponds to E(A = 0), II(A — 
1), A(A = 2), <3>(A = 3) etc. symmetries. There are several methods to solve equation (9), but we 
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adopt an integral-equation technique by converting the differential equation (9) into an integral 
equation 33 , 

Plk M = Mkr)S " „ + XI ^(( r > 7,1 ) ( rl io (r')dr' , (10) 

v Jo 

where the Green’s function is defined by 

Ge(r,r') - k~ 1 [r) ( (kr)j t (kr') - jt(kr)rj t (kr')], (11) 

in which je(kr) and rji(kr) are Riccati-Bessel functions. Note that in the expansion ( Eq. 8 ) of 
polarization term, we have only A = 0 and 2 terms. 

The CO molecule is a polar molecule which needs special attention in a BF adiabatic-nuclei 
theory where the forward DCS and a t are undefined 35 . In this respect, we employ the multipole- 
extracted-adiabatic-nuclei (MEAN) scheme of Norcross and Padial 36 , in which the DCS for the 
J —> J 1 rotational transition are given as 

£[<W; oo)! ! p | 2 (Bh. - b f J*)p,(cosi>), (12) 

It X — 0 




where the first term is the usual closed form for the (JJ’) rotational excitation DCS in the space- 
fixed first-Born-approximation (FBA); kj and kj> are respectively the wavevectors in the initial and 
the final channels; C (...) is a Clebsch-Gordan coefficient; l t is the angular momentum transferred 
during the collision; B\, i t are the DCS expansion coefficients and B^f t A are the corresponding 
quantities in the FBA evaluated in the BF frame of reference. The channel vectors are related by 
the relation, 

k 2 j-k 2 j, =B[J'(J' + 1)~ J(J + l)), (13) 

where B is the rotational constant of the CO molecule. Finally, the expressions for total ( ) 

and the momentum transfer ( ) cross sections are evaluated from equation (11) for any (J J') 

transition. Total ( summed over all final rotational states J f ) integrated ( a t ) and momentum 
transfer ( cr m ) cross sections can easily be obtained from 

E j jt 

a t or m 1 ( 14 ) 

J 1 

However, for a proper comparison with experiment, we average the a t and <r m over the Boltzmann 
distribution of rotational states at 300 K ( represented as (a t ) and (<r m ) ). This is quite easy since 
in the present energy region the sum over in equaton (14) is insensitive to J. 

Finally, for a polyatomic molecule, we employ totally a different set of computer codes to 
obtain scattering parameters. The details are given elsewhere^ 13 > 37 \ For a polyatomic target, it is 
again convenient to make use of the ANA and set up the scattering equations in the single-center- 
expansion scheme under the close-coupling formalism. The equation for the scattered positron 
function P(r) is the same as given in Eq. (7), however, now the single-center-expansion scheme 
is quite different. The angular basis functions belong to the irreducible representation (IR) of the 
molecular point group. The P( r), V st and V po i are expanded around the center-of-mass (COM) of 
the molecule in question. For example, the P( r) is expanded as 13 , 

Plrl^r-SSr’WlWlf), (15) 

thpfi 
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where the are symmetry-adapted angular functions belonging to a particular irreducible rep- 
resentation pfi of the molecular point group. The static potential V 3t is given by 

N M 

Vjt(r) = /|tM 2 l r - ^\~ 1 dr 1 dr 2 • • ■ ^ Z t \r - K\\ 3 , ( 16 ) 

^ j= 1 i= 1 

where I s the target ground state wavefunction given as a single Slater determinant of one- 
electron spin orbitals <£ a (r), N is the total number of electrons and M the number of nuclei in the 
molecule. 

Finally, the scattering amplitude for with initial direction k and the final direction f is 

m.t) = § - *«•*“•). (17) 

where the sum being over l, h,£',h',p and p. The S matrix is related with the K matrix in the usual 
fashion. The amplitude ( Eq. 17 ) is defined in the BF frame of reference. In order to transform 
it into the space-fixed (SF) or the laboratory frame coordinate system, we employ the standard 
technique in terms of rotation matrices 13 , 3?(a,/3,7) ( the (a/?7) are the three Euler’s angles ). If 
r' represents the coordinates of the positron with respect to SF coordinate system, the transformed 
amplitude, /( k.r'), is employed to determine the rotationally inelastic transition amplitude under 
the ANA theory, i.e., 

/(i-f) = (i|/(k.r';a/? 7 )|f), 

where {i| and (f| are respectively the initial and final rotational eigenfunctions. The total elastic 
cross sections are obtained by summing over all final rotational states and averaging over all initial 
states. The expressions for the differential, integral and momentum transfer cross sections are given 
for general non-linear polyatomic molecules in Ref. 13 . 

In the present CH4 case, we do not face any convergence problem in the summation over 
various angular momentum quantum numbers. In fact, even the DCS’s can be obtained easily with 
proper convergence within a reasonable size of the scattering matrix. For more details about the 
actual numerical parameters we recommend our previous paper 29 . In the present results on the 
positron-CH 4 collisions, we have kept the same single-center expansion and K-matrix parameters 
as described in Ref. 29 . 

Jain and Thompson 29 used three different approximations for the polarization interaction; all 
the three models were exactly the same as employed for electron scattering^ 38 -40 F However, the 
most successful was the one based on the second— order perturbation theory under the Pople Shofield 
method 41 in which the distortion in each molecular orbital is the same. The non-adiabatic effects 
were included via the non-penetration criterion of Temkin 42 . This electron polarization potential 40 
( to be denoted by JT ) has so far been quite successful in positron-CH 4 collisions. Unfortunately, 
the JT potential has never been employed for any atomic system; however, it has been employed for 
the e + -N 2 case 43 with fair success. The other two EPP used by Jain and Thompson 29 were based 
on the asymptotic form multiplied by the cut-off function. Although the use of non— penetrating 
scheme for the positron case may be questionable, however, we argue here that the non-adiabatic 
effects in e + case are not much effective due to relatively smaller local kinetic energy of the impinging 
positron in the vicinity of the target. Therefore, we emphasize here that in any polarised-orbital 
variational approach, the non-adiabatic correction may not be taken very seriously. In the findings 
of Elza et al 8 , these non-adiabatic effects seem to change the results; however, it is hard to draw 
any conclusion since they introduce fitting parameters which makes physics less clear. 


77 



Fig. L Various polarization potentials for the e + -Ar system. The po- 
tentials PCOPl , PCUP2 and ECOP are shown , respectively , by solid, 
dashed and long-dashed lines. The various notations are explained 
in the text. In the inset are shown the total optical potentials ( sum 
of repulsive static and attractive polarization potentials ) with vari- 
ous models: solid line , V 5t + V P £ 0P2 ; dashed line , V Jt + V r ^ OPl ; 

iong-dasbed iine, V 5t + V p ^ op . 


IV. Results and Discussion 


First, we display the new PCOP terms in Fig. 1 for the e^-Ar system. We have considered 
two forms of the PCOP model; one, the correlation potential defined by Eq. ( 6 ) ( to be denoted by 
PC0P2 ) and two, the correlation energy itself ( Eq. 3, to be denoted by PCOPl ). Both the PCOPl 
and PCOP 2 terms are plotted in Fig. 1 . Also shown in this figure is the corresponding ECOP 
potential. We see a significant difference between the ECOP and PCOP curves. In general, the 
PCOP is stronger then the ECOP approximation. This simply means that the e + — e“ correlation 
energy is stronger, thus giving rise to a more attractive polarization potential. It seems realistic 
since the is expected to distort the target charge cloud deeper due to strong correlation of 
electron and positron particles. A similar situation exists for molecular targets ( CO and CH 4 ) ( 
not shown ). 
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Fig. 2 Differential cross sections for the e + -Ar elastic collision at 8.5 eV. 
Present theory: solid curve, PCOP2 model; dashed curve, PCOPl 
model; long-dashed curve, ECOP model. The experimental points 
are taken from Refs. 44 and 45. 


The total optical potential shown in the inset of Fig. 1 is a sum of repulsive static and attractive 
polarization terms; thus there is a zero-potential point and an attractive well. The position of the 
zero-point potential and the shape of the attractive well decide the penetration and sign of various 
partial waves. Here the role of polarization interaction is important and the low-energy scattering 
is strongly influenced by these cancellation effects not present in the electron scattering case. Thus 
the form of the polarization potential in the zero-potential and the attractive well region is very 
crucial to determine the scattering process. 

Fig. 2 illustrates the DCS for the e + -Ar system at 8.5 eV alongwith the ECOP results and the 
measurements of Refs. 44-45. We see a significant qualitative improvement represented by these 
PCOP curves. The dip in the experimental DCS around 50° is neatly reproduced by the new model, 
while the ECOP dip occurs at smaller angle (30°). We have not shown other calculations 46-49 due 
to their semi-empirical nature. We have seen similar agreement between theory and experiment 
at other energies ( lower and higher than 8.5 eV ) also 16 . In order to further see the success of 
the PCOP model at lower energies, we have calculated the scattering length in the zero energy 
limit. The value of the PCOP scattering length is —4.89 (au), which compares very well with the 
experimental value 50 of —4.4 ± 0.5 (au); the ECOP model gives this value only 1.7 (au). It is thus 
quite clear that in this low energy limit a true positron polarization interaction makes big difference 
in the scattering parameters. 
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Fig. 3. Positron-CO ( <r t ) cross sections in various theoretical models ( 
including the R-matrix results 3 , shown as dash-dot-dot curve ) and 
experimental observations ( circles, Kwan et a l 51 ; triangle, Sueoka 
and Mori 52 ). Our PCOP results are shown using experimental value 
of the dipole moment ( sohd curve ) and also the theoretical value of 
the dipole moment ( dashed curve ). The ECOP data 16 are plotted 
as the dotted curve. 


We now discuss our e + -CO calculations on the (cr t ) parameter. Fig. 3 shows the present PCOP 
(<T() alongwith the ECOP 28 , R-matrix 4 and the experimental^ 51 -52 ^ results. We have shown two 
versions of the present POOF ( Eq. 6 ) model: one, with the present theoretical value of the dipole 
moment (D =0.099 a.uT), shown as dashed curve in Fig. 3 and two, by using the experimenta l 
value (D=0.044 au), shown as a solid line. Below 4 eV, the use of theoretical value of the dipole 
moment makes large changes in the { <r t ). The R-matrix results do not compare well with the 
measu reme n ts; o ne reason being that they employ their theoretical di pole momen t value in the 
MEAN approximation and the second reason is related with their polarization force which still 
needs to be improved 53 , The inclusion of the polarization force makes Targe changes ( about a 
factor of two or more ) in the pure static results. Ii seems that the ECOP approximation is better 
at higher ( E > 2 eV ) energies; however, this conclusion may not be true as the two sets of DCS’s 
differ significantly and we expect that the PCOP DCS’s are better than the cooresponding ECOP 
cross sections 1 ' . 
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Fig. 4. Positron-CO DCS at 2 eV in the PCOP ( solid curve ) and ECOP 
( dashed curve ) 26 polarization models. 


In order to see the qualitative differences between the ECOP and PCOP models for the positron 
collisions with CO molecules, we have shown in Fig. 4, the differential cross sections at 2 e\ . We 
see a qualitative difference between the two curves of Fig. 4 . The dip in the ECOP approximation 
occurs at lower angles while the PCOP dip in the DCS’s curve occurs at somewhat larger angles; 
the difference between the positions of the two dips is about 20°. The disagreement between these 
two polarization model is seen at all angles ( Fig. 4 ). It is interesting to note that the total cross 
sections at this energy is almost same in both the PCOP and ECOP approximations. Thus, the 
integral cross sections are sometimes confusing and therefore a theoretical model should be judged 
from the angular functions which are more sensitive to model potential results. It would be very 
interesting to see the position of the dip in an experimental investigation. So far, we believe that 
the PCOP dip is more realistic. 
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Fig. 5. The positron-CH 4 <x t cross sections using the PCOP , ECOP 28 
and JT 29 theoretical models. All the curves axe labelled with respect 
to all these approximations. The experimental data are taken from 
Refs. 54 and 55. 


From Fig. 3, we thus see that the PCOP model is quite successful in reproducing the ex- 
perimental data, particularly at low energies where the ECOP calculations totally fail. At higher 
energies ( above 3 eV ) where the ECOP approximation seems to work well, the difference betw een 
the two sets ofTDCS is significant see Fig. 4 ]; for example, at 2 eV, the positions of the dips in 
the ECOP and PCOP curves occur at 40° and 50° angles respectively. . ....... . 

Bight now there are no other DOS’s ( theoretical or experimental ) available for comparison 
for the positron— CO elastic collisions. It would have been interesting to have a comparison between 
the PCOP and the R-matrix angular functions. It is possible that the DCS’s for the positron-CO 
system be measured in the laboratory in future. Only then one can conclude finally about the 
usefulness of the PCOP model in the e + -CO case. In passing, we would like to mention that 
at further lower energies ( below 1 eV ) the difference between various DCS calculations may be 
dramatic due to a stronger dependence of the collision dynamics on the polarization/correlation 
effects. - 
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Fig. 5. The positron-CH 4 <r t cross sections using the PCOP . ECOP 
and JT 29 theoretical models. All the curves are labeUed with respect 
to all these approximations. The experimental data are taken from 
Refs. 54 and 55. 


From Fie. 3, we thus see that the PCOP model is quite successful in reproducing the ex- 
perimental data, particularly at low energies where the ECOP calculations totally fail. At higher 
energies ( above 3 eV ) where the ECOP approximation seems to work well, the difference between 
the two sets of DCS is significant ( see Fig. 4 ); for example, at 2 eV, the positions of the dips in 
the ECOP and PCOP curves occur at 40° and 50° angles respectively. 

Right now there are no other DCS’s ( theoretical or experimental ) available for comparison 
for the positron-CO elastic collisions. It would have been interesting to have a comparison between 
the PCOP and the R-matrix angular functions. It is possible that the DCS s for the positron CO 
system be measured in the laboratory in future. Only then one can conclude finally about the 
usefulness of the PCOP model in the e--CO case. In passing, we would like to mention that 
at further lower energies ( below 1 eV ) the difference between various DCS calculations may be 
dramatic due to a stronger dependence of the collision dynamics on the polanzation/correlation 

effects. 


83 




Angle (Deg) 


Fig. 6. DOS's for the e + -Cfl 4 collisions at 6 eV using the POOP, ECOP 
and JT 29 models for the polarization potential. The points ( + ) are 
the static only results. 


Finally, in Fig. 5, we demonstrate our e + -CH 4 <r t values along with experimental points. Also 
shown in this Fig. 4 are the calculations using EPP models of Jain and Thompson 29 and Gianturco 
et al 28 . The measured data are taken from Dababneh et al 54 and Sueoka and Mori 55 . Again we 
see that the new PCOP gives very encouraging results as compared to other theoretical models 
based on the EPP. However, if we look into the corresponding DCS ( Fig. 6 ) at 6 eV, the three 
models give totally different type of results. We need experimental data on the DCS jn order to 
see which polarization approximation describes the collision properly. It is clear from Fig. 6 that 
even if there is some agreement in the total cross sections, the DCS’s may be quite different both 
in quality and quantity. 


V. Conclusions 

We conclude that a true positron polarization approximation is essential to investigate the 
low-energy positron collisions with atoms and molecules. Even if the EPP and PCOP total cross 
sections are identical, the DCS differ significantly. We are in the process to test the sensitivity of 
the Z e ff parameter with respect to EPP and PCOP polarization models. We mention here that 
in the positron scattering, the question of nonadiabaticity may not be as serious as in the case of 
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electron collisions: the reason being that the positron being repelled by the nucleus of the target is 
not speeded up in the vicinity of the target; thus reducing the non-adiabatic effects. 

In this article, we have discussed a new parameter-free positron polarization potential which 
is obtained by short range e + — e~ correlation energy and the correct long range behaviour of the 
polarization potential. It should be realized that the inclusion of correlation/polarization effects in 
positron collisions is a very hard problem. In model approaches, such as the one discussed here, 
one has to compromise with the local and energy-dependent form of the polarization against its 
non-local and non-adiabatic nature. The use of model polarization potentials will still continue 
due to their simplicity and significant success even at the DCS level. We emphasize that all the 
positron polarization models should be tested with respect to differential cross sections and the 
very low energy parameters such as the scattering length and the annihilation paremeter. 
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